{p, TT^) correlations in central heavy-ion collisions at 1 -^ 2 AGeV 

A.B. Larionov ^, W. Gassing, M. Effenberger, U. Mosel 
Institut fiir Theoretische Physik, Universitdt Giessen, D-35392 Giessen, Germany 

Q^ -. Abstract 

The proton - charged pion correlated emission is studied in the reactions Au 

o ; 

^ ■ (1.06 AGeV) + Au, Ni (1.06 and 1.93 AGeV) + Ni and Ni (1.97 AGeV) + Cu 

^^ ■ within the BUU approach. The associated invariant mass distributions are 
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shifted to smaller energies with respect to the free A (1232) mass distribution 
due to kinematical reasons. We find that the existing and partly conflicting 
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r^ • experimental data do not allow to draw definite conclusions on the in-medium 

o 



modification of the A(1232). 



I. INTRODUCTION 



The properties of the baryon resonances in hot and dense nuclear matter are the subject 
'k>< I of recent theoretical [|l]-P and experimental 0-§] studies. The most important question is, 

p^ \ how the centroid energies and the widths of the resonances are modified in nuclear matter. 

In Ref. the invariant mass spectra of correlated {p, vr^) pairs were measured in colli- 
sions of Ni+Cu at 1.97 AGeV. It is demonstrated in |0, that the peaks of the spectra are 
shifted to smaller invariant masses by about -50 MeV with respect to the free A mass; this 
shift increases with the collision centrality. Furthermore, in Refs. ||^J^ the mass distributions 
of the baryon resonances excited in Au+Au and Ni+Ni central collisions at energies between 
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1 and 2 AGeV were deduced experimentally on the basis of two methods: (i) by defolding 
the pt spectra of charged pions, and (ii) from invariant masses of {p, n^) pairs. The peaks 
of the A-mass distributions extracted in Ref. ^ are shifted on average by —83 MeV for 
Au+Au and —63 MeV for Ni+Ni collisions with respect to the free A-resonanceQ. 

In this work we have performed BUU transport calculations (see Refs. |Q,^ for a general 
review and Ref. 0] for the description of the latest version of the BUU model) for collisions 



of Au+Au at 1.06 AGeV, Ni+Ni at 1.06 and 1.93 AGeV and Ni+Cu at 1.97 AGeV in order 
to study the invariant mass spectrum of correlated {p, vr^) pairs produced by resonance 
decay. The aim of this work is twofold: (i) to study the mechanism of the correlated pair 
emission, and (ii) by comparison with the experimental data try to get information on the 
in-medium modification of the A-resonance. 

The structure of this paper is as follows: In Sect. II the space-time picture of the corre- 
lated {p, vr^) pair emission is studied. In Sect. Ill the final correlated {p, vr^) invariant mass 
spectra are presented and compared to thermodynamic model predictions. The possibility 
to extract these spectra within the mixed event technique |TT| is demonstrated and a com- 



parison to the experimental data of Refs. |§-§[ is performed. A summary and conclusions 
are given in Sect. IV. 

II. SPACE-TIME PICTURE OF THE CORRELATED (j9, TT^) PAIR EMISSION 



All calculations performed in this work involve the BUU model of Ref. |T^, including a 



momentum-dependent mean field with a compressibility modulus K = 280 MeV as well as 
vacuum spectral functions for all resonances. We shall call this parameterset as "standard" 



^ The mass shifts reported in Ref. Q are given relative to the free A peak of 1.210 GeV, which 
arises from the ;B-function of Ref. [^ (see also dashed line in Fig. 14). In our work, the mass 
distribution of the free A resonance is peaked at 1.224 GeV (see solid line hne in Fig. 14) and we 
discuss all mass shifts with respect to this value. 



below for brevity. Modifications of spectral functions will be discussed explicitly. Here we 



only briefly describe the most important physical inputs of the BUU model (see Ref. |]T0 
for details). 

The mechanism of pion production implemented in the BUU model includes two steps 
(besides a small contribution from direct A^A^ -^ NNtt processes): In the first step the 
A (1232) and higher baryon resonances are excited in inelastic nucleon-nucleon collisions: 
A^A^ — i> NR, where A^ stands for a nucleon and R for a resonance. At a collision energy l-h2 
AGeV dominantly the A^A^ -^ A^A(1232) channel determines the inelastic part of the total 
nucleon-nucleon cross section. The masses of the produced resonances are chosen according 
to their vacuum spectral functions 

A(M) = ^ ^^^^^) (1) 

^ ' vr (M2 - M2„iJ2 + M2r2(M) ' ^ ^ 

where Mpoie is the pole mass and r(M) is the total mass- dependent decay width. In the 
particular case of the A-resonance Mpoie = M^ = 1.232 GeV, while 

r(M)^r.(^)'^|±4 (2) 

\qrj M (3^ + q^ 

is the free space mass- dependent A-width with Fa = 0.118 GeV, [3r = 0.2 GeV; q is the 
pion momentum in the rest frame of the A and g^ is the value of g at M = Ma. 

In the second step the resonance decays into a pion and a nucleon: R -^ ttN. The 
resonance decay at a given time step is simulated by Monte-Carlo using the free space mass- 
dependent decay width r(M). The Pauli blocking for the outgoing nucleon is taken into 



account in our calculations (c.f. Ref. [0). The life time of the A-resonance at its pole 
mass with respect to the pion emission is, therefore, ta^ttAt = l/FA^TrAf ~ 2 fm/c, where 
^A~>nN is the in-medium decay width of the A-resonance, which is less than r(M) due to 
Pauli blocking ||T^. Furthermore, the produced pion can be reabsorbed: vrA^ -^ A. Since 
the average life time of the A is much less than the characteristic time (~ 40 fm/c) of a 
central heavy-ion collision at 1^2 AGeV and the pion mean-free-path is quite small {X-^ ~ 1 
fm), a chain of processes Ai -^ 7!'iNi, ttiN2 -^ A2, A2 — > t^sN^, ... develops. Moreover, the 



pion and nucleoli rescattering on nucleons strongly reduces the probability for the decayed 
A to be "visible" . Therefore, only a very small part of A-resonances excited in the nucleus- 
nucleus collision, which is not absorbed in the process AA^ -^ NN, can be observed by 
looking at their decay products - i.e. correlated (vr, A^) pairs. We define a proton-pion pair 
as being correlated, if this pair originates from the same resonance and both the proton and 
the neutron don't rescatter anymore. In other words, we will always consider only observable 
correlated pairs. For the collisions under study we have obtained numerically on average l-h4 
(p, vr") and 2-^6 (p, vr^) correlated pairs per event in a calculation without acceptance cut 
(Table 1). After filtering through the central drift chamber (CDC) acceptance of the FOPI 
Collaboration (for Au+Au and Ni+Ni reactions) the number of correlated pairs decreases 
by a factor of 3-i-4. The CDC acceptance (see Refs. |]7|,p^) was simulated by selecting only 
particles in the interval 32" < Qiab < 150° and, for 7r+, additionally the cut piab < 0.65 
GeV/c was applied. In Table 1 we present the numbers N^°^^--. and A^fi'^'V) of correlated 
(p, TT") and (p, vr+) pairs and corresponding ratios r"^ = N?p^^±J NT"-^±y Here A'L^^i) is the 
number of "random" pairs, i.e. those pairs which are composed of a proton and a pion from 
the same event, but not emitted from the same resonance: 

KZ^) = N,K^ - Ar(— ±) , (3) 

where Np and A^^± are the total numbers of emitted protons and pions, respectively. The 
ratios r^ are quite small (~ 10"'^ -ir 10^^) which makes the separation of the true signal from 
the background (see Sect. Ill B) difficult. 

In Fig. 1 we show the time evolution of the nucleon density profile and of the radial 
distribution of the correlated {p, vr^) pairs emitted during the time intervals from tj — 1 fm/c 
to tj + 1 fm/c, where ti=14, 18, 28 and 38 fm/c, for the Au+Au system at b=0. Most of the 
pairs are emitted from the periphery of the system, where the density is low (cf. Fig. 3). 
We see from Fig. 1, that the intensity of the emission reaches a maximum at t ~ 20 fm/c 
and then slowly decreases. This is consistent with the pion production rate, which reaches 
a maximum at t=20 fm/c (cf. Fig. 11 from Ref. [|l^). Thus, pairs are emitted dominantly 



at 20 < t < 30 fm/c. 

In Fig. 2 the radial dependence of the temperature is presented for a central collision 
of Au+Au at various times accompanied by the radial distribution of the correlated {p, vr^) 
pairs. The local temperature was determined on a cubic lattice of 1 fm grid size introduced 
in the center-of-mass system of the colliding nuclei. A finite temperature Fermi distribution 
has been adjusted such as to obtain the correct baryon density and average baryon kinetic 
energy in the local rest frame of the matter element in the lattice cell. During the time 
period of intensive pair emission, 20 < t < 30 fm/c, the temperature at the periphery 
of the system decreases from 70 MeV to 30 MeV. This implies, that the description of 
the freeze-out of correlated {p, ii'^) pairs by some universal value of the temperature is quite 
schematic. The same conclusion can be obtained with respect to the freeze-out density, which 
is varying strongly during the period of intensive pair emission (Fig. 1). Nevertheless, one 
can estimate some average value of the freeze-out density and temperature for a comparison 
with macroscopic thermodynamic models (cf. Refs. [p|,p^). Fig. 3 shows the distributions 
of the correlated {p, it) (tt = tt^, tt^, 7r+) pairs in the freeze-out temperature-density plane for 
different reactions. It is interesting that these distributions are very similar for all colliding 
systems regardless of the system size and, moreover, the collision energy. We always see 
a doubly peaked structure: at Tf^. ~ and pf,o. ~ 0, that corresponds to the decay of 
A-resonances in vacuum, and at Tfo. ~ 23 MeV and pfo. ~ 0.14po, where po = 0.16 fm~^ 
is the equilibrium nuclear matter density. In Table 2 the mean values and dispersions of 
the freeze-out parameters are given for all reactions. The distributions are quite broad[| 
and the mean values of the freeze-out density and temperature are essentially higher than 
the maxima of the actual distributions in Fig. 3: {pf.o.) — 0.41po and {Tf,o.) — 39 MeV. 
The baryon chemical potential extracted within the thermal model (see next section) at 



^ The numerical tables with the distributions of the {p, n) pairs in the freeze-out temperature- 
density are available from the authors as data files. 



Pb = 0.41po and T = 39 MeV is hb = 889 MeV. These freeze-out parameters are different 
from the thermal model analysis of Cleymans et al. for Au + Au collisions at 1 AGeV (Ref. 
0): i.e. Tf,o. = 50 MeV, ^b = 850 MeV. 



III. INVARIANT MASS DISTRIBUTIONS OF CORRELATED {p, TT^) PAIRS 

First we discuss the results of the direct analysis of the correlated {p, vr^) pairs which 
can be unambiguously identified in the BUU calculation. In Fig. 4 the invariant mass 
spectrum of correlated {p, vr^) pairs is shown (solid line) for Au+Au central collisions. This 
spectrum was extracted by stopping the time evolution at t = 40 fm/c. All the correlated 
pairs present in the system at this time plus the pairs produced by the forced decay of 
resonances were counted. We have checked that stopping the time evolution at later times 
does not influence the extracted spectrum of the pairs within the accuracy of our statistics, 
since for Au + Au central collision at 1.06 AGeV only 2% of all baryon-baryon collisions 
happen after 40 fm/c. In order to understand, how various time intervals are contributing 
to the total invariant mass spectrum, we show also in Fig. 4 the partial spectra given by 
the pairs emitted at t < 20 fm/c and at t < 30 fm/c. We see that, in agreement with the 
discussion in the previous section, the total spectrum is dominantly composed from pairs 
emitted in the interval 20 < t < 30 fm/c and that the low-mass pairs come primarily from 
very late times. 

A. Analysis within the thermodynamical model 

We have analysed the final spectra of correlated {p, vr^) pairs using a simple thermody- 
namical model including nucleons, A(1232) resonances and free pions. For a given tempera- 
ture T and baryon density pb the baryon chemical potential ps can be exctracted from the 
equation: 

Pb = Pn + Pa , (4) 



where pn and pA are the densities of nucleons and deltas: 

- A [ ^^P ^ (K\ 

P""' J ^2^ exp(ep-/iB)/T + l ' ^ ' 

p. = 16 / .M^(M,/^-^-^^-^ (6) 

with ep = \/M'^ + p^. In Eq. (^ ^(M) is the spectral function of the A-resonance given 
by Eq. (|1|). The mass distribution of the A-resonance at finite chemical potential fiB and 
temperature T is 

'"''^ A{M)IW I ,fl„ , \ ,^ , . (7) 



dM "■ ' J (27r^)3 exp (cp - Pb)IT + 1 

where V is the freeze-out volume. 

In Fig. 5 we show the invariant mass (p, vr) distributions from the BUU calculations 
and the A-mass distribution Eq. (|^) at various temperatures for the baryon density ps = 
{Pf.o.) = 0.41po- For comparison, the spectral function A{M) (^ is also shown (dot-dashed 
lines). The high- mass tail of the spectral function (M > 1.3 GeV) is populated very weakly. 
Furthermore, the resulting (p, vr) invariant mass spectra have a shape quite different from 
the function A{M) due to the presence of the integral over the Fermi distribution in Eq. (^. 
The peak of the distribution dN^/dM is slightly shifting to higher masses with increasing 
temperatures. We see that the BUU invariant mass spectrum is close to the thermal one for 
T ~ 40 ^ 50 MeV. 

A better understanding of the BUU results is achieved if we weight the thermal model 
calculation with the distribution function 9^A^pair(7/.o.,P/.o.)/<9T/.o.9p/.o. of the (p, vr) pairs 
(Fig. 3). Then the invariant mass spectrum of the pairs reads as follows (see Appendix for 
a derivation): 

it J J dTdp y ^miv(M'\§Ea 

dM 



-^^ = dT dp ^J^'Z 00 — • (8) 

dM J J '^ dTdp J dM'T{M')^^^MM:^LEl 



Mn+M^ 

The distribution (^ is shown by solid lines in Fig. 5. In Eq. (^) we calculated the inte- 
grals for limits T = 5 -7- 100 MeV, p = 0.1 -r- 2po. In this way the pairs emitted early at 



nonequilibriuin (high density) and late (low density) stages of the heavy-ion collision were 
removed from the analysis. The weighted distribution fits the BUU spectrum better than 
a calculation with fixed temperature and density using Eq. (|^). Note that Eq. (^ has no 
free parameters: it only uses as an input the distribution of the emitted pairs at freeze-out 
temperature and density. Therefore, the mass distribution of the emitted pairs from central 
collisions is consistently described with a local thermal equilibrium assumption in agreement 
with the analysis in Refs. [0,||. 

B. Extraction of the correlated (p, vr^) pairs by background subtraction 

The experimental extraction method of the correlated {p, vr^) pairs is based on back- 
ground subtraction from the event-by-event spectrum [|6|-[^,pl]] . A standard way of the back- 
ground construction is the event mixing technique [|lT|. We have, therefore, also prepared 



mixed BUU events taking protons and pions from different events. In the case of Au+Au 
collisions all events were taken at b=0 fm. For the reaction Ni+Ni we have taken events 
for b = 1,2 and 3 fm, and a mixed event can be composed of two events at different impact 
parameters^. We have checked that an additional restriction of equal impact parameters in 
both events does not influence the result in this case, since the collision dynamics changes 
strongly only for b~ 5 fm. For the reaction Ni-|-Cu a mixing of b = 1,2. . .5 fm events 
was performed selecting event pairs with equal impact parameters only. Then the difference 
spectrum agrees with the spectrum of real pairs. However, for the Ni+Cu reaction, the 
mixing of the events without imposing any condition on the impact parameters produces a 
difference spectrum peaked at about 1.15 GeV, i.e. at much smaller invariant mass, since 
the background spectrum gets shifted to higher invariant masses. The conclusion is that 
the difference spectrum resembles the spectrum of real pairs only if both events selected 



^The background construction was done only for the beam energy of 1.93 AGeV in the case of 
Ni+Ni collisions. 
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for mixing have impact parameters within 3 fm. Note, that in our BUU study the reaction 
plane is fixed and, in distinction to the analysis of the experimental data in ||^, we did not 
rotate our events around the beam axis. 

In Figs. 6-8 the total event-by-event and background spectra (a,c) and their differences 
(b,d) are shown for the systems Au+Au, Ni+Ni and Ni-|-Cu. For the event mixing procedure 
we have prepared sets of 30000, 66000 and 24000 BUU-events for Au+Au, Ni+Ni and 
Ni+Cu reactions, respectively. In order to reduce the computational time, these events 
were calculated with the Coulomb interaction switched off (for a discussion of the Coulomb 
effects see the next section). The number of mixed events is 10 times more for each reaction. 
The real correlated pairs, selected event-by-event, were extracted in parallel from the BUU- 
events. 

The difference spectra (points with errorbars in Figs. 6-8 (b,d)) reveal a clear correlation 
signal at the invariant mass M ~ 1.2 GeV, in agreement with the spectra of the correlated 
pairs (solid lines). 

The general agreement between the spectra of real {p, vr^) pairs and the difference spectra 
opens the possibility to identify the spectrum of real pairs given by BUU with the experi- 
mental difference spectrum. In this way one can get rid of large statistical errors, which are 
mostly due to the background constructed from mixed BUU events. At the same time, when 
doing such a comparison, we neglect some systematical deviations of the difference spectrum, 
which depends on the adopted background construction procedure, from the spectrum of 
real correlated pairs, which is unambiguous. For instance, we see from Figs. 6-8 that, in 
distinction to the real spectrum, the difference spectrum can show a negative correlation 
and an increased high invariant mass tail. Therefore, we will concentrate only on the gross 
structure (like the peak position and the width) of the calculated and measured spectra. 



C. Comparison of BUU with the experimental data 

In Figs. 9-11 we show the invariant mass distributions of the correlated {p, vr") - (a) and 
(p,7r+) - (b) pahs for Au+Au at 1.06 AGeV (Fig. 9), Ni+Ni at 1.06 AGeV (Fig. 10) and 
Ni+Ni at 1.93 AGeV (Fig. 11) reactions in comparison to the experimental data of Ref. 
1^. The experimental acceptance of the CDC was taken into account in our calculations. 
The theoretical curves are averaged over impact parameter in the range 6 < 3 fni, that ap- 
proximately corresponds to the PM5 multiplicity bin (see Ref. [0). Results with a standard 
parameterset of the BUU model are shown by dotted lines. These results are obtained taking 
into account the Coulomb interactions between charged particles in the BUU calculation. 
However, it was supposed in the standard calculation, that the momenta of all particles 
are frozen after 40 fm/c and, therefore, the residual Coulomb energy was neglected. The 
effect of the residual Coulomb energy was, furthermore, taken into account by rescaling the 
momenta of particles as: 

p ^ Kp, K= {2mUcoui/p^ + 1)'/' , (9) 

where Ucoui is the Coulomb energy of a particle at t = 40 fm/c. The residual Coulomb 
energy shifts the spectra of {p, vr") and {p, vr"*") pairs to smaller and larger invariant masses, 
respectively, by about 5-i-lO MeV (dashed lines). 

The complete role of the Coulomb interactions can be seen now, e.g. from a comparison 
of the solid line in Fig. 6d with the dashed line in Fig. 9a. We see that the Coulomb effects 
are shifting the peak of the (p, vr^) spectrum in the Au+Au central collisions by about 
-50 MeV. The corresponding shift of the {p, tt'^) spectrum is about +20 MeV for the same 
reaction. At higher beam energies of about 2 AGeV the role of the Coulomb interactions 
becomes negligible for the Ni+Ni and Ni+Cu systems. 

We have also performed a calculation considering the width of the A-resonance selfcon- 
sistently, i.e. taking into account the Pauli blocking of the decay A — » A^vr, the absorption 
AA^ —>■ NN and the rescattering AA^ —>■ AN. This dynamical width of the A-resonance 
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then has been used in Eq. (|T|) for the spectral function (see Ref. ^7^ for details). Moreover, 
since quite high central densities ~ 2.5po are reached in central heavy-ion collisions at 1 ^ 2 
AGeV, a medium modification of the matrix element for the process AA^ -^ NN simulating 



the effect of a direct three-body absorption (see also Ref. [O) was included as follows 



\Man^nn\' = f 1 + 3^ j \MTn^j,^\' ■ (10) 

The solid lines in Figs. 9-11 show the calculations including the selfconsistent A- width (wsc) 



and the density-dependent matrix element (ddme) of Eq. (|ToD- It gives essentially broader 
invariant mass spectra and shifts the peaks to smaller invariant masses, however, still not 
enough to explain the FOPI data on {p, 7t~^) pairs. 

One observes, nevertheless, a rather good overall agreement of the wsc+ddme calculation 
to the FOPI data on (p, vr"). For the Ni+Ni coUisions at 1.93 AGeV (Fig. 11a), the {p, vr") 
data reveal a double-humped structure with an additional peak at low invariant mass, which 
can be explained by the decay A -^ pn^ (c.f. Ref. 0). This decay creates a narrow peak 
at the invariant mass of 1.116 GeV. We obtained on average ~ 0.1 (p, vr^) pairs per event 
in central collisions of Ni+Ni at 1.93 AGeV and of Ni+Cu at 1.97 AGeV (see below) due to 
the A decays. The channels A^A^ -^ AKN and A^A^ -^ T.^KN followed by TP -^ A7 were 
taken into account in the calculation of the A production. The A contribution to the (p, vr") 
spectrum is shown by the dot-dashed line on top of the wsc+ddme curve in Fig. 11a. This 
contribution was parameterized by a gaussian of width a = 6.4 MeV/c^ determined by the 
bin size 20 MeV/c^ of the calculated spectrum. 

Figs. 12,13 show the results for the Ni+Cu reaction at 1.97 AGeV in comparison to the 
data from Ref. 0. For this reaction we have performed only the wsc+ddme calculation. 
The calculated results are not filtered through the EOS-TPC acceptance. Again, we can 
observe, however, a good agreement with the (p, vr") data taking into account the A decay. 
We note, that the relative A contribution is larger in the Ni+Cu reaction than in the very 
similar Ni+Ni reaction. This is caused by the 5+6 times reduction of the number of {p, vr^) 
pairs due to the A decays after filtering through the CDC acceptance, while the number 

11 



of {p, vr^) pairs produced by the A decays gets reduced only by a factor of 3-i-4 (see Table 
2). We attribute this different reduction to a smaller directed A flow than proton flow (c.f. 



Ref. |T8[)- For (Pj^t^) pairs (Figs. 12b and 13), the calculations still overpredict the peak 
position of the A++ by about 10 MeV, but the overall agreement with the EOS-TPC data 
on {p, 7T~^) pairs is better than for the FOPI data. 

It was shown in Ref. p, that the peak of the {p,tt'^) invariant mass distribution shifts 
to lower invariant masses with the centrality of the collision. Fig. 13 shows, that also in 
the BUU calculations the same effect is present since in peripheral nucleus-nucleus collisions 
the emitted proton-pion pairs are mostly due to decays of A-resonances excited in energetic 
first-chance nucleon-nucleon collisions, where the thermal picture discussed in Sect. Ill A 
does not apply. We should remark that our statistics is rather poor for peripheral collisions. 
Therefore, the discussed shift in our calculations in Fig. 13 is more relevant for the average 
value of the distribution than for the peak position. 

D. Comparison of the BUU-|-thermal calculations to the experimental data 

We have studied the predictions of the thermal model using as an input the distribution 
of (p, vr) pairs at the freeze-out temperature and density produced by BUU (see Eq. (|D). 
This hybrid approach offers an easy possibility to see the influence of the spectral function 
in the thermal part of the calculation on the observed A-mass spectrum while retaining in 
BUU the bare spectral function A of Eq. (|l]). Besides the calculations with the bare spectral 
function A in the thermal part, we have performed a thermal model analysis replacing A 
in Eqs.(P),(0) by the derivative of the vrA^-scattering phase shift in the P33 channel with 
respect to the center-of-mass energy of the pion and nucleon (Ref. 0): 

B{E,^) = 2^^^^. (11) 



The i3-function ( [TT| ) can be interpreted as a level density of the irN system (cf. Ref. \T^\ and 
Refs. therein). The thermal model employing the i3-function (^) gives a better agreement 
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with the experimental vr^ spectrum for the Au+Au reaction at 1.06 AGeV as shown in Ref. 
^. To clarify the reason we compare in Fig. 14 the spectral function A (solid line) with 
the weight function B/2n (dashed line) (see also Fig. 2 in Ref. 0). The i3-function yields 
an increased contribution at smaller invariant masses. 

Fig. 15 demonstrates the results of the hybrid BUU+thermal calculations based on Eq. 
(1) for various systems in comparison to the data from Refs. PJ^. For this comparison, we 
have selected the data on {p, vr"'"), since: (i) the statistics for {p, tt^) is always better than for 
{p, vr"), and (ii) the {p, vr^) signal from the A^"^ is not contaminated by the decays of higher 
baryon resonances and by the A decay. The solid lines in Fig. 15 show the BUU + thermal 
calculation with the ^-function. This calculation strongly overpredicts the peak positions 
in the case of Au+Au and Ni+Ni reactions (see also dotted lines obtained within the BUU 
only in Figs. 9-11 b). The agreement is somewhat better in the case of the Ni+Cu reaction. 
For the i3-function (dashed lines), the peak of the theoretical spectrum of the pairs shifts 
to smaller invariant masses while the width increases. Thus, the agreement of the thermal 
calculation with the FOPI data is improved by the i3-function. The residual discrepancies 
can be further diminished by accounting for the Coulomb effects in the calculations (c.f. Ref. 
1^). However, the S-function leads to a worse agreement for the Ni+Cu reaction, which can 
be explained reasonably well only within the BUU (wsc+ddme) calculation (dotted lines in 
Fig. 15). This calculation, however, is still not consistent with the FOPI data on (p, tt"*"). 

IV. SUMMARY AND CONCLUSIONS 

In this work a study of correlated {p, vr^) pair emission from central heavy-ion collisions at 
energies of l-i-2 AGeV has been performed within the BUU transport model. In agreement 
with the data |^ less than 1% of the total number of the {p, vr^) pairs are correlated. 
Our calculations give ~ 25% of emitted pions in correlations with protons (the rest of 
pions are produced either directly or their correlations are destroyed by proton and/or 
pion rescattering) for Au+Au central collisions at 1.06 AGeV. This value is lower than the 
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one reported in Ref. of > 50% of pions correlated with protons. Since the total pion 
multiplicity is also overpredicted by the BUU model for the Au+Au reaction by about a 
factor of 1.7, the number of correlated pairs per event turns out to be close to the data 
again. 

The correlated pairs originate from A-decays in low-density regions during an intermedi- 
ate stage (t = 20 -T- 30 fm/c) of the collision. The calculated invariant mass spectra of these 
pairs have a thermal shape; however, the real (kinetic) temperature in the system changes 
quite strongly during the period of emission. This result is related to a weak sensitivity of 
the shape of the invariant mass spectra to the temperature in the region T = 40 -h 55 MeV, 
which follows from a simple thermodynamical calculation (see text and Fig. 5). We would 
like to stress, that the freeze- out of the pairs happens during an extended period of time, com- 
parable with the time scale of the heavy-ion collision itself. Therefore, the thermodynamical 
state of the system is changing quite strongly during this period. 

The BUU calculation with the selfconsistent treatment of the A- width and the medium- 
modified matrix element for the process A^A -^ NN (see Sect. Ill C) gives a good overall 
agreement with both the FOPI and EOS-TPC data on the {p, vr") invariant mass distribu- 
tions. The calculated {p, vr^) spectra are shifted to higher invariant masses with respect to 
the FOPI data, but are in reasonable agreement with the EOS-TPC data. It is, therefore, an 
open question if some additional effects influencing the A propagation and decay in hot and 
dense nuclear matter are needed to account for the remaining shift in future dynamical cal- 
culations. In particular, the approach taking into account the vacuum vrA^ scattering phase 
shift (Refs. p|,|20[|) gives a better agreement with the FOPI data, as demonstrated within 
the thermodynamical model in Ref. [0 and in Sect. Ill D of the present work; however, it 
gives a worse description of the EOS data from the BEVALAC This conflicting situation 
calls for new experiments on {p, vr) correlations. 
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APPENDIX 

In this Appendix we derive Eq. (§) for the mass distribution of emitted pairs. 
The density of emitted pairs in the space invariant mass - freeze-out temperature - 
freeze-out density is (for brevity we drop lower indices at Tfo. and Pf.o.)'- 



dMOTdp ' 'J dMdTdp 

where — dMdTd ^^ ^^^ density of A-resonances in the same space as a function of time, 
r(M) is the decay width of Eq. (j^), and W is the surviving probability of an emitted pair 
assumed to be independent on invariant mass, temperature and density. We have: 

^%4M= / mwvmJ^t^^^MM^^, (13) 

dTdp J ^ 'J dMdTdp ^ ' 

MN+M-n 

Assuming local thermal equilibrium we can write: 

d^NA{M,T,p,t) _ d^NA{T,p,t)dNA{M,T,p) 
dMdTdp ~ dTdp dM 

where 



(14) 



dNA{M,T,p) _ dNA{M,T,p) 
dM dM 



oo 



,,,,SN.(M',T,p) 

dM' 



-1 



(15) 



is the mass distribution of A-resonances (see Eq. (|^) normalized to 1, which depends only 
on the local temperature and density. 

Substituting Eqs. ([T^),(p!5|) into Eqs. (|T2|) and (0), we obtain the expressions: 
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dMdTdp - ^^^^' ^)^r(M) — , (16) 

%f^^a(P,T) 7 ,MWVm'-^^^^ (17) 



with 



Using Eqs. (p!^) , (|17|) and the relation 



dM J J '^ dMdTdp 
it is straightforward to obtain Eq. (| 
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Table 1 Numbers of correlated {p, vr^) pairs N?°^^±^ and ratios of numbers of correlated and 
random (p, vr^) pairs r^^'> (see text for definitions) for central collisions. In the case 
of Au+Au and Ni+Ni reactions, the two numbers are given for N?°'^^±^ separated by a 
comma: before and after filtering through the CDC acceptance. The calculated values 
of r*^^) are practically insensitive to the filtering. The experimental values of r*^^) from 
Ref. are given in brackets. For Ni+Cu collisions all results are unfiltered. 



System 


Energy 


TVfcorr 


r(-) 


]\Tcorr 


r(+) 




(AGeV) 




(%) 




(%) 


Au+Au 


1.06 


4.0, 1.2 


0.1 (0.6+0.2) 


5.7, 2.1 


0.3 (0.75+0.25) 


Ni+Ni 


1.06 


1.2, 0.4 


0.4 (0.75+0.25) 


2.4, 0.7 


0.9 (1.0+0.3) 


Ni+Ni 


1.93 


2.5*, 0.6 


0.4 (0.6+0.2) 


4.9, 1.2 


0.9 (1.05+0.3) 


Ni+Cu 


1.97 


2.7* 


0.3 


4.9 


0.7 



including 0.1 (p, vr ) pair due to A decays 



Table 2 Mean values and dispersions of the freeze-out density and temperature for central 
collisions. 



System 


Energy 


(Pf.o.) 


CTp 


{Tf.o.) 


ax 




(AGeV) 


(Po) 


(po) 


(MeV) 


(MeV) 


Au+Au 


1.06 


0.41 


0.39 


39.0 


21.0 


Ni+Ni 


1.06 


0.45 


0.42 


39.4 


25.6 


Ni+Ni 


1.93 


0.40 


0.42 


37.8 


24.9 


Ni+Cu 


1.97 


0.41 


0.42 


37.9 


24.9 



19 



FIGURE CAPTIONS 

Fig. 1 Radial distributions of the correlated {p, vr^) pairs emitted during the time interval 
At = 2 fm/c (histograms) and the baryon density profiles (solid lines) at various times 
for a central collision of Au+Au at 1.06 AGeV. Pairs are selected without acceptance 
cuts. 

Fig. 2 The same as Fig. 1, but with the kinetic temperature shown by solid lines. 

Fig. 3 Freeze-out temperature-density distribution d'^NpairiTj o., Pf.o) IdTf o.dpf.o. of the 
correlated (p, vr) pairs for central collisions of various systems. The distribution (in 
a.u.) at a given point {Tfo., Pf.o.) is proportional to the size of the box. No acceptance 
cuts were applied. 

Fig. 4 Invariant mass distribution of the correlated (p, vr^) pairs emitted at the time inter- 
vals t = O-i-20 fm/c (dotted line), t = 0^30 fm/c (dashed line) and the total spectrum 
after induced decay of residual A's at t = 40 fm/c (solid line) for a central collision of 
Au+Au at 1.06 AGeV. No acceptance cuts were applied. 

Fig. 5 {p, tt) (tt = 7r~, n^, 7r+) invariant mass spectra from the BUU calculations for central 
Au+Au collisions at 1.06 AGeV for all pairs (filled circles) and for pairs emitted at 
Tf.o. = 5 + 100 MeV and pf,o. = 0.1 ^ 2po (filled squares) in comparison to the thermal 
model calculation at the temperature T = 39 MeV and the baryon density pb = 0.41po 
(dotted line), T = 55 MeV and pb = OAlpQ (dashed line), and weighted with the 
distribution function d'^NpairiTf o.i Pf.o)/dTf o.dpf,o. according to Eq. (||) (solid lines). 
The free spectral function A{M) (Eq. (|l])) is shown by the dash-dotted line. Curves 
for fixed temperatures and the free spectral function are normalized to the number of 
pairs emitted at Tf o. = 5 ^ 100 MeV and pf,o. = 0.1 + 2po- No acceptance cuts were 
applied in selecting the pairs. 

Fig. 6 The left panels show the spectra of all (p, 7r+) - (a) and (p, vr^) - (c) pairs extracted 
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event-by-event (solid line) and by event mixing (dashed line) from the BUU calculation 
for Au+Au at 1.06 AGeV. The right panels present the corresponding differences 
between event-by-event and mixed-event spectra of (p, vr"*") - (b) and (p, tt") - (d) 
pairs (points with errorbars connected by dashed line) and the spectra of real pairs 
from the same BUU calculation (solid line). The particle selection was performed 
including the CDC acceptance. 

Fig. 7 The same as Fig. 6, but for Ni+Ni collisions at 1.93 AGeV. 

Fig. 8 The same as Fig. 6, but for Ni+Cu collisions at 1.97 AGeV. Particles are selected 
in full 4:71 acceptance. The A-decay contribution to the (p, vr") spectra is dropped. 

Fig. 9 Invariant mass spectra of (p, tt^) (a) and (p, vr^) (b) pairs from central Au+Au 
collisions at 1.06 AGeV in comparison to the data from Ref. 0. Dotted lines - 
standard BUU parameterset, dashed lines - standard BUU with Coulomb final state 
interaction, solid lines - BUU calculation with selfconsistent width of the A-resonance 
and a density-dependent matrix element (see text for details) plus Coulomb final state 
interaction. All spectra are normalized to unity. 

Fig. 10 The same as Fig. 9, but for central collisions of Ni+Ni at 1.06 AGeV. 

Fig. 11 The same as Fig. 9, but for central collisions of Ni+Ni at 1.93 AGeV. The dot- 
dashed line shows the A decay contribution to the spectrum of (p, vr") pairs added to 
the A° decay spectrum calculated with selfconsistent width of the A-resonance and 
density-dependent matrix element. 

Fig. 12 The same as Fig. 9, but for Ni+Cu collisions at 1.97 AGeV. The BUU calculation 
is performed with selfconsistent A-width and density-dependent matrix element for 
the impact parameter region b=0 -^ 10 fni. The final state Coulomb interaction is 
taken into account. The A decay contribution to the spectrum of {p, vr") pairs added 
to the A° decay spectrum is shown by the dot-dashed line. The experimental data are 
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from Ref. 0. 

Fig. 13 Invariant mass spectra of (p, vr"*") pairs for Ni+Cu collisions at 1.97 AGeV at dif- 
ferent centrality. The BUU calculations (performed with selfconsistent A-width and 
density-dependent matrix element) are presented by histograms for b=2,4,6,8 and 10 
fm from the uppermost to lowermost panel. The Breit-Wigner fits to the data from 
Ref. [^ are shown by dashed lines for the multiplicities M > 55, 45 < M < 55, 
35 < M < 45, 25 < M < 35 and M < 25. The spectra are normalized to unity. 

Fig. 14 The spectral function A of the free A (solid line) and the weight function B/27r 
(dashed line) versus invariant mass. 

Fig. 15 Thermal invariant mass distributions using functions A (solid lines) and B (dashed 
lines) weighted according to Eq. (§) in comparison to the data on {p, vr"*") pairs from 
Ref. (Au+Au and Ni+Ni central collisions) and from Ref. (Ni+Cu collisions 
with multiplicity M > 55). The dotted lines show the BUU (wsc+ddme) results for 
central collisions. All curves are normalized to unity. 
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